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Abstract 

The recently observed IceCube PeV events could be due to heavy dark matter (DM) decay. In 
this paper, we propose a simple DM model with extra U{l)x gauge symmetry and bridge it with 
standard model particles through heavy right-handed neutrino. The Dirac fermion DM y with 
mass ~ 5 PeV can dominantly decay into a dark Higgs (b), the SM Higgs (/i) and a neutrino 
(n). If the lifetime of x is ~ 0(10^®) sec, the resulting neutrino flux can fit data consistently. 
The neutrino flux from x <fhv in our model is softer than the one predicted from x n/i, for 
example. We also discuss a possible mechanism to produce DM with the right relic abundance. 
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I. INTRODUCTION 


Recently, the IceCube Collaboration has reported the detection of 37 neutrino events 
with energy between 30 TeV — 2 PeV, of which three have energy above 1 PeV [1-3]. 
According to the recent analyses [4-7], the three-year data are consistent with equal fluxes 
of all three neutrino flavors and with isotropic arrival directions. However, the neutrino flux 
required to £t the data in lOOTeV — PeV range is around 10“® GeV cni“^s“^sr“^ per flavor, 
and rejects a purely atmospheric explanation at 5.7a. Therefore astrophysical and/or new 
physics explanations have been pursued for the origin of these high energy neutrinos. 

Possible astrophysical sources are involved with supernova remnants (SNR) [8-12], ac¬ 
tive galactic nuclei (AGN) [13-16], and gamma-ray bursts (GRB) [17, 18], all of which 
assume some specihc emission spectra due to different production environments. In a model 
independent analysis, the best-£t power law spectrum from the IceCube analysis [1] is 
Eld^u/dEy ~ 1.5 X 10“^(i7i,/100 TeV)“°'^ cm“^ s“^sr“h It is not very straightforward to fit 
such a spectrum by astrophysical sources. In all cases, extragalactic sources are needed due 
to the isotropic feature and galactic constraints [19-22]. Then identifying such astrophysical 
sources will be crucial for further understanding of IceCube events. 

Dark matter (DM) and other new physics interpretations have been also investigated 
in various ways or models. Heavy DM might decay into SM particles that give energetic 
PeV neutrinos [23-37] or it could decay into some light DM particles which interact with 
nucleons and mimic neutrino events [39, 40]. The resulting neutrino flux would still be 
consistent with isotropy so far, since galactic and extragalactic DM contribute at the similar 
order. One unique feature of DM explanation is that there should be sharp energy cut-off in 
the neutrino spectrum, which could be tested by future data. Also, a possible gap around 
400 TeV ~ 1 PeV, although not statistically signihcant yet, motivated considerations of new 
interactions, two-component flux and leptophilic DM decay [32, 41-47]. 

In this paper, we propose a simple DM model to explain the IceCube PeV events. A 
dark sector with new U{l)x gauge symmetry is introduced and can have connection with 
the standard model (SM) sector through neutrino-portal interactions as well as the Higgs 
portal interaction. Fermionic DM (y) has ~ PeV mass and mostly decays into three-body 
hnal state with dark Higgs (</>), SM Higgs (h) and active neutrino. The produced neutrinos 
from primary y decay and the secondary h and 0 decays can explain the observed PeV event 
spectra, while the atmospheric and astrophysical neutrinos are included for the low-energy 
part. 

This paper is organized as follows. In Section. H we introduce our DM model with dark 
U{l)x gauge symmetry, heavy right-handed neutrino portal and Higgs portal interactions. 
In Section. HI we outline the general formalism for calculating neutrino flux from galactic and 
extragalactic DM decay. In Section. IV we present both total and differential decay width 
for the relevant three-body decay in our model and compare the numerical results with 
IceCube data. In Section. V we discuss a possible mechanism to generate the correct relic 
density for DM and direct/indirect detection constraints. Finally, we give our conclusions. 


If involving of dark halo substructure, Ref. [38] showed annihilating DM scenario may also be possible. 
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II. MODEL 


We consider a dark sector with a dark Higgs field d> and a Dirac fermion DM x associated 
U{l)x gauge symmetry. Their U{l)x charges are assigned as follows 


Qx) — ( 1 ’ !)• 

We begin with the following renormalizable and gange invariant Lagrangian inclnding jnst 
one singlet right-handed (RH) nentrino N and one lepton flavor (more Ns and/or flavors 
can be easily generalized): 

C =£sm + ]^Ni$N - {UnxN^N + ylHN + h.c.^ - ^ sin 

+ - R(<h, H) + x X - [fx^X + h.c.), (2.1) 

where L = {u l)'^ is a left-handed (LH) SM SU{2) lepton donblet, H is the SM Higgs 
donblet, = 9^X,^ — d^X^ is the field strength for U{l)x gange field X^, Fy^ is for SM 
hypercharge f/(l)y, and e is the kinetic mixing parameter. Two types of Ynkawa conplings, 
y and /, can be taken as real parameters, ignoring CP violation for simplicity. We define 
covariant derivative as = df^ — igxX^. Since we are interested in explaining the IceCnbe 
PeV events in terms of DM x decay, we shall take ~ PeV. Other parameters in onr 
model are free variables. 

The scalar potential V of this model is given by 

X = A„ - y) + A*f, (t'-J - y) + A* (^4>*4> - , (2.2) 


Both electroweak and dark gange symmetries are spontaneonsly broken by the nonzero 
vacnnm expectations valnes of H and <h: {H) = (O, vh/V^) , (*h) = v^/\/2 . Here vh — 
246GeV is the same as SM valne bnt might be taken as a free parameter. In the nnitarity 
gange, we can replace the scalar fields with 


A A 

y/2\VH + h[x) 


and <h —)■ 


+ (p{x) 

V2 ' 


(2.3) 


Note that h and (j) shall mix with each other thanks to the Higgs-portal operator, (the X^h 
term) Throngh this mixing, </> can decay into SM particles. Another important mixing 
happens among three nentral gange bosons, photon A^, and X^. Snch a mixtnre wonld 
enable an extra mass eigenstate Z/ (mostly X^) to decay SM fermion pairs. Then DM 
X scattering off nnclens is possible by the Z' exchange, and the cross section essentially 
depends on e, mz'- It is easy to choose small e, or heavy masses to evade the constraints 
from DM direct detection [54]. 

When the right-handed nentrino N is mnch heavier than x, we can integrate it ont and 
obtain an effective operator, 

^X^H^L + h.c., (2.4) 

rriN 


^ A similar setup with different dark charge assignments has been considered for the AMS02 positron 
excess [48]. One may also use discrete symmetries, see Ref. [49] for example. 

^ The Xh 4 , term can also help to stabilize the electroweak vacuum [50-53]. 
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which would make x decay possible but long lived. After spontaneous gauge symmetry 
breaking, we have several higher dimensional effective operators from the aforementioned 
operator Eq. (2.4) as follows: 


V4>Vh _ 
- 

rriN 


v<t> 

rriN 


xhjy, 


Vh 

niN 




1 

niN 


X(j)hu, 


(2.5) 


with the common factor 


yl 

2 


for all these operators. If kinematically allowed, all the above 


operators induce y decays into different channels with fixed relative branching ratios. Within 
the heavy y limit, the mass operator yz/ in Eq. (2.5) would 

induce a tiny mixing between y and u with the mixing angle (3 approximately given by 


[3 ~ 


yf v^vh 

2 m^rrij^ 


( 2 . 6 ) 


Then the gauge interactions for y and p will generate the decay channels. 




(2.7) 


with their branching ratios being proportional to ^ vjj : : 2v‘^. Two dim-4 operators, 

xhp and X(j)P, would lead y to the following decays. 


X ^ hp, (j)P, (2.8) 

with their branching ratios being proportional to ~ vl : vj,. It is also straightforward to get 
the following relation for the branching ratios, 

Br{x —t (j)p) : Br{x —t Z'p) ~ 1 : 1 . (2.9) 


Therefore, all the decay branching ratios are basically calculable and completely fixed in 
this model Note that the decay modes with Z' or 0 are unique features of DM models 
with dark gauge symmetries 

Another interesting phenomenon in this model is that three body decay channel x (php 
is dominant over all other channels when S> Vfp: 


Ta (x -t (t>hp) _1 

Ts (x ^ hp, (t)p) ~ 167r2 vl + v% 


( 2 . 10 ) 


since we actually have an enhancement from heavy even though there is a phase space 
suppression from three-body final states. There are another three-body decay channels that 
are equally important: 

X —t <\)jZ' h p, <PIz' -\- Z -\- p^ 4^1 z' -\- -j-, 


with branching ratios 1:1:2 due to the Goldstone boson equivalence theorem. In the 
following, if not otherwise stated explicitly, we use x 4hp to represent all these channels 
and in numerical calculations we take all of them into account. 


^ This is also true in the model for the AMS02 positron excess [48] . 

® This is also true of three-body decays of DM discussed in the following paragraph. 
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To give a rough impression for the relevant parameter ranges, we can perform an order- 
of-magnitude estimation (complete formulas and details of calculation are given in the Ap¬ 
pendix) : 


r3 (x -» 4'hv) ~ 




967r^ 

y^rriN J lO^^sec 

~10-^®GeV-\ 

TJIn 



( 2 . 11 ) 

( 2 . 12 ) 


which shows the required value for the combination yf/rriM. If one additionally assumes that 
N should be responsible for neutrino mass through the usual Type-I see-saw mechanism, 
then we have y 10 ^y/uiTv/PeV. Therefore we would have y 1 and / ~ 10 for 
rriN ^ 10^"^GeV and, y ~ 10“^ and / ~ 10“^® for tun PeV. In any case, / is very tiny and 
seems unnaturally small. However, it is still natural a la ’t Hooft [55] since taking / = 0 
enhances symmetries of the theory, namely DM number conservation In later discussion, 
we shall take y, f,mN as free parameters, unless specihed. 


III. NEUTRINO FLUX FROM DM DECAY 


The neutrino flux from dark matter decay is composed of galactic and extragalactic 
contributions which are equally important as we shall see below. Galactic neutrino flux at 
kinetic energy E from DM decay in our Milky Way dark halo is given by 




dE„ 


E^=E 



dK 

dEy 



(3.1) 


where P* is partial width for decay channel i, dNl/dE^ is the neutrino spectrum at pro¬ 
duction, r' = -f — 2rQr cos 9, r is the distance to earth from the DM decay point, 
Tq ~ 8.5kpc for the solar system and 9 is the observation angle between the line-of-sight and 
the center of the Milky Way. For the galactic DM density distribution, we use the following 
standard NFW prohle [56], 


p° ('•') = p@ 


r®i 

1 -h Tq/Tc 

. r'. 

1 -|- r '/tc 


(3.2) 


with parameters Tc — 20 kpc and pg ~ 0.4 GeV/cm^. For decaying dark matter the flux 
is not very sensitive to DM density prohle, so our discussions and results will still apply if 
another different prohle is used. 

We can also get the extragalactic or cosmic contribution from a formula similar to the 
above one, by taking cosmic expansion into account, namely the red-shift ehect [24]: 




dE„ 


E„=E 


Pc^X 

47rm. 


Er. 


dz dN^ 


n dE„ 


(3.3) 


El/ — (l-\-z)E 


where E' is red-shifted to as = (1 -1- z)E, the red-shift 2 ; is dehned as 1 4- ^ = ao/a 
with present scale factor oq being normalized to 1, the critical energy density pc = 5.5 x 


® The DM current = xi^X is conserved in the limit / —>■ 0. 
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FIG. 1. Neutrino spectra from DM x decay with ~ 5 PeV and lifetime = l/P ~ 2 x lO^^s. 
The left panel shows individual contribution of different final states from y’s decay, ly (blue dot- 
dashed curve) and h/(f) (red dashed curve), respectively. The right panel presents the galactic (blue 
dashed curve) and extragalactic (red dot-dashed curve) neutrino flux. 


10 ®GeV/cm^ and ~ 0.27 is DM x’s fraction. The Hubble parameter "H is related to its 
present value through 


77 — T/o'x/Da + Dm(l -\- z)^ -|- f2j.(l -|- z)'^, 

Da, Dm and D^ are energy fractions of dark energy, all matter, and radiations, respectively. 
We shall use the latest results from Planck [57] for numerical evaluation. 


IV. NUMERICAL RESULTS 


To compute the neutrino flux from DM decay, we hrst need to calculate the total and 
differential three-body decay width for x ^ (p + h + u. In the heavy % limit, we have obtained 
the total width 


T ~ ^ 


yf 

76871^ KrriN 


and normalized differential decay widths 

1 dr 


T dEi, 

1 dr 


1 dr 


TdE. 


~ 2AEf^/rrf^, Q < E^ < m^/2, 

~ l2Ei, - Ej,) /mj, 0 < E,, < mJ2, 

^ I 2 E 4 . {m^ - lm\, 0 < < m^/2. 


(4.1) 


(4.2) 

(4.3) 

(4.4) 


The details of the calculation are given in the Appendix where complete formulas with 
nonzero mass parameters are also presented. The above differential widths are essential 
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ingredients to get the final neutrino fiux. For example, z/s from different decay final states 
are given by 


dN 

~dE 


{x —>■ u) 


1 dr dN, {E,) 
r dEx dE 


dEx, 


(4.5) 


where x = i>,h, W, Z, Z', 0. Note that dNi, {Ex) /dE in the integrand can be calculated with 
Pythia [58] or PPPC4DMID [59]. 

In Fig. 1 we show the neutrino spectra. Just for illustration, we choose the mass of DM 
X around 5 PeV and its lifetime 2 x 10^®s. The neutrino spectra are multiplied by E/ in 
order to account for the energy dependence of the neutrino-nucleus cross section so that it 
might be compared with IceCube data more easily. The blue dot-dashed curve indicates 
the spectrum from final v in y’s decay and red dashed one marks h/0’s contribution. For 
simplicity, the mass of 0 has been chosen to be just as the SM Higgs mass. However, other 
choice does not affect our result much, since in the left panel of Fig. 1 we see that in the 
high energy part h/ 0 ’s contributions are basically negligible and it is the high energy part 
that explains the IceCube PeV events. In the right panel, we include the red-shifted effects 
and show both galactic (blue dashed) and extragalactic (red dot-dashed) contributions. 

Next, we compare our model predictions with IceCube three-year data [1]. To parame¬ 
terize the possible astrophysical neutrino fluxes at low energy, we consider either a broken 
power law (BPL) or unbroken power law (UPL) [32], 


2 dd>bkg 

" dE, 

2d<Fbkg 


Jr 


BPL 


E, 


lOOTeV 


-71 


exp 


J 


UPL 


VlOOTeVy 


-72 



(4.6) 

(4.7) 


where the first one has an exponential cut-off at energy scale Eq which is chosen to be 125 
TeV in agreement with the SNR results [ 8 ] . 

We illustrate in Fig. 2 with two different astrophysical fiux choices. The low energy data 
are best fitted by varying Jq and the spectral index 7 , but the high energy PeV data points are 
fit by DM decay. In the left panel, we have used = 4.1 x 10“®GeV/cm^/s/sr and 71 = 
0(red dot-dashed curve). In the right panel, we have used = 1.3 x 10“®GeV/cm^/s/sr 
and 72 = 0.7(red dot-dashed curve). DM’s contributions and total fiux are labeled with 
purple dashed and blue solid curves, respectively. As we can see in the figure, our model can 
agree with the PeV data. Also a gap could appear around 400 TeV although the current 
data can not tell existence of the gap is statistically significant. The feature for the DM 
decay spectrum is that there would be a sudden drop around E, ~ '^x/2- As more data are 
accumulated, it should be possible to test our model in the future. 

We should note that in our discussion Jq and 7 are just adjusted visually to be consistent 
with the low-energy data points. Dedicated investigation would require global fitting, which 
is beyond our scope here. Just for comparison, IceGube [60] gives the best-fit parameters of 
a single unbroken power law for neutrinos energies between 25TeV and 2.8PeV without DM 
contribution, Jq = x 10“®GeV/cm^/s/sr and 7 = 0.5 ± 0.09. 

In Fig. 3, we also compare our model with the preliminary updated results [61] based 
on IceGube 4-year data which has already filled the gap a bit. Here we have only shifted 
to J®^'" = 5.6 X 10“®GeV/cm^/s/sr and ~ 1.5 x 10^®s (left panel), and J(^^^ = 2.1 x 
10“®GeV/cm^/s/sr and ^ 2 x 10^®s (right panel). In Fig. 4, we illustrate two cases, one 
with = 8 PeV, ~ 1.7 x 10^®s, and the other = lOPeV, ~ 1.5 x 10^®s. 
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FIG. 2. Neutrino flux from DM x’s decay with ~ 5PeV and lifetime = l/F ~ 2 x 10^®s 
and IceCube Data [1]. The left (right) panel used a broken (unbroken) power law (BPL) for 
astrophysical neutrino flux with a red dot-dashed curve. DM’s contributions and total flux are 
labeled with purple dashed and blue solid curves, respectively. See details in the text. 
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FIG. 3. Same as Fig. 2 but with preliminary updated results based on 4-year data [61], with 
= 5.6 X 10“®GeV/cm^/s/sr and ~ 1.5 x 10^®s (left) , and = 2.1 x 10“®GeV/cm^/s/sr 
and ~ 2 X 10^®s (right). 

There are some crucial differences between our model and some others in the literature. 
For example, the authors in Ref. [23, 29] considered the effective operator, yLHx with 
y ~ 10“^°, which induces mainly two-body decay of DM y, 

X ^ z/h, 

In this scenario, the neutrino spectrum shows that there should be no gap between 400 TeV ~ 
1 PeV [26]. Our model predicts that the dominant decay mode are 

X —y (f)/z' -\- h -\- 4^1 z' -\- z -\- 4 /T + ^^5 
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FIG. 4. Same as the right panel of Fig. 3 but with = 8PeV, r,^. ~ 1.7 x 10^®s (left), and 
= lOPeV,~ 1.5 x 10^®s (right). 


which is a consequence of U{l)x dark gauge symmetry and the dark charge assignments of 
the dark Higgs and dark matter fermion y. The neutrino spectra from primary y decay and 
the secondary decays of h and 0 have different shapes and could account for the possible 
gap. However, we should note that the current data can not favor one over another yet due 
to its low statistics. Also the neutrino flux in our model is softer than the one predicted in 
Ref. [23, 29], for example. 

In Ref. [32], leptophilic three-body decay induced by dimension-six L^lpL^x was consid¬ 
ered with global U{1) or flavor symmetries. Besides the neutrino spectrum difference, 
our model involves an additional gauge boson which mediates the DM-nucleon scattering, 
and could be tested by DM direct searches. 

Our scenario is also different from those in which DM decay is also responsible for the 
low-energy flux [24]. The DM lifetime in Ref. [24] should be around 2 x lO^^s, as mainly 
determined by the low energy part of events. This is partly due to the reason that the 
branching ratio into neutrinos and bb there should be about 10% and 90%, respectively, to 
account for the possible gap. On the other hand, in our scenario 1/2 of the decay channels 
have prompt neutrinos. Another main difference is that three-body-decay usually gives 
broader spectra at PeV range than two-body-decay considered in Ref. [24], but more data 
is required in order to discriminate this difference. 

Assuming the dark photon is much heavier than Z, the DM-nucleon scattering cross 
section can be roughly estimated as 




rrir 


m 


X 


sin^ e X 10 ^®cm^. 


(4.8) 


10 ^®cm^ is the typical cross section value for SM Z-mediating DM-nucleon process. Com¬ 
paring it with the direct detection bound for lOOGeV DM, we should have 


< 10 ^^cm^ X 


m. 


lOOGeV’ 


(4.9) 
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for heavy ~ 5 PeV. This can be easily satished, for example, with mx ~ TeV and 
sin e < 0.1. 

V. RELIC ABUNDANCE AND CONSTRAINTS 

In onr above investigation, we have not discnssed the relic abnndance for DM y yet. Since 
y is very heavy, the nnitary bonnd on its annihilation makes the DM y impossible to be 
thermally prodnced and a non-thermal process is needed (see Ref. [62] for a recent review 
on snch topics). Here, we discnss one possible non-thermal prodnction mechanism for DM 
y in onr model. We assnme that the dark Higgs (j) ‘ once shared a common temperatnre 
with SM particles and had a thermal distribntion when its temperatnre T was larger than 
m^. Here we do not specify the mechanism how 0 reached snch a temperatnre; it conld be 
dne to reheating after inflation or some heavy particle decays. 

In the thermal bath, y conld be prodnced throngh 0 -f- 0 —>■ y -f y, whose thermal cross 
section is given by 



Here we considered the mx 3> T case only. We can calcnlate the y’s yield, = n^+^/s, 
where s ~ g^s{T)T^ is the total entropy density in the Universe ( g^siT) ~ 100), 



(5.1) 


In the above derivation we have nsed the Hnbble parameter "H ~ ^/gJ^T‘^/g*{T) ~ 
100 is the total effective nnmber of degree of freedom when 0’s temperatnre is T and Planck 
mass Mpi = ~ 4.2 x 10^® GeV. Since the yield has a positive power dependence 

on temperatnre, y is mostly prodnced at 0’s highest temperatnre 


max’ 



(5.2) 


Reqniring y give the correct relic density, we have a relation 



(5.3) 


which pnts a constraint on / and m^v. 



(5.4) 


When mx > Tmax > '^x ~ PeV, we are able to give a lower bonnd for / 

I/I > 10-^ 


^ It should be <i> precisely since at high temperature symmetries are not yet broken, but it will not affect 
our discussion. 
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One more thing we can infer from Eq. (2.11) and (5.4) is that in this production mecha¬ 
nism y would be too small such that this right-hand neutrino N can not be fully responsible 
for active neutrino mass and mixing angles. This is because, in order to explain the active 
neutrino mass, we should have 


rriN 


~ O.leV ^ ~ 10 ® 


mN \ 


(5.5) 


which can not be satisfied simultaneously with Eq. (2.11) and (5.4). This is not a problem 
for this model since we can expect there are additional right-handed neutrinos Ni and they 
can couple to LH with large yi but to with tiny /*, so that they are just responsible for 
active neutrino mass and mixing angle but not for DM y’s production. 

DM direct detection can constrain DM-nucleon scattering cross section, whose value in 
our model is determined by the kinetic mixing parameter e, gauge coupling gx and the mass 
of Z'. For GeV-TeV DM, there is already plenty of viable parameter space to evade such a 
constraint, see Ref. [63] for example. For PeV DM, the constraint is even relaxed due to the 
low number density, see Eq. 4.9. Indirect detection from positron, anti-proton and 7 rays 
also constrain DM y’s decay lifetime, see Refs. [64-66] and references therein for example. 
We have checked that ~ 10^® sec is still allowed by all such constraints. 

As an illustration, in Fig. 5 we show the expected gamma-ray flux from DM decay with 
~ 5PeV and lifetime ~ 2 x 10^®s. We have included the prompt gamma-rays and 
those from inverse compton scattering (ICS) of charged particles on CMB, starlight and 
dust-rescattered light. For extragalatic contribution, we have taken absorption factor into 
account. As we can easily see, the flux is well below the current constraints from Fermi- 
LAT [67] and KASCADE [ 68 ]. Our results are also consistent with the recent investigations 
about gamma-ray constraints on the lifetime of PeV DM in different scenarios [36, 69], 
> 3 X lO^T 


VI. CONCLUSION 

In this paper, we have proposed a dark matter (DM) model that can explain the IceCube 
PeV events in terms of DM decay. The model is based on an extra U{l)x dark gauge 
symmetry which is spontaneously broken by a dark Higgs held, <h. One crucial bridge 
between DM y and standard model (SM) particles is established by heavy right-handed 
(RH) neutrino portal interactions. This heavy neutrino can induce DM decays into SM 
particles, including the light neutrinos. 

The dominant decay channel of DM is the three-body hnal state with SM Higgs, dark 
Higgs, and neutrino (y —)■ (f) + h + p) (and other channels due to the Goldstone boson 
equivalence theorem), and not the usual two-body decays such as y —)■ Zz/, etc.. This 

is a unique feature of the present model based on U{l)x dark gauge symmetry and the RH 
neutrino portal interactions. We have calculated both total and differential decay width 
to evaluate the galactic and extragalactic neutrino fluxes. We have found that neutrino 
flux from these decay products can agree well with the IceGube spectrum. Together with an 
astrophysical flux for lower energy events, we are able to £t IceGube data around 100 TeV ~ 
2 PeV if we assume DM mass is about ~ 5PeV and its lifetime is ~ 2 X 10^8 sec. 
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FIG. 5. The gamma-ray flux from DM decay with ~ 5PeV and lifetime ~ 2 x lO^^s, 
confronted with constraints from Fermi-LAT [67] and KASCADE [68] data. 
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APPENDIX 


Here we show the complete differential decay width for x ^ h + (p + u. Throughout the 
calculation, we work in the rest frame of y, so y’s momentum is {m^, 0,0, 0). For unpolarized 
X, we have 


dV 


1 1 
(27r)^ 8m^ 


J2\M\^dE,dEf,= 

pol 


1 



E^dEjydEh^ 


( 6 . 1 ) 


where we have used the averaged, squared matrix element. 






2 

m^Ey. 
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Then we get 


dr 


E,, 


dEy (27r) X^rriNJ Je 


yf \ 


2 ^Smax 


dEh 


E,. 


(27r) \2mN 




2m. 


2E,. 


yf 


{2n)m^ \2m]\f 

yf 


Ef-ml^Ef-ml 


(27r) \2mN 


(with E^ < m^/2). 


(6.2) 

(6,3) 


In the last line we have used the heavy limit, 3> m/i, 3> mpu and ~ 0. Some 
dehnitions are listed below, 


Ek = 


+ m\ — 


(piy 


2m^ 


^Ib = i.Pa +Pbf , rUab = \l{Pa+Pbf, 
and some other kinematic variables. 


+ ^ - mj^ + ml 

2mh4> ’ '' 

=^l + ml- 2m^E^, 

where E'^ and E^ are the energies of 0 and u in the mhip rest frame, respectively. From the 
above differential decay width, we can easily get the total width 


I ral / yf V 

24 (27r)^ \2mN J 


(6.4) 


The differential decay width as function of Eh or E^ can also be calculated similarly. For 
example. 


dF 

dEh 


yf Y 


(27r)^ \2mN J J e^' 

1 (ji\ 


EydEy 


'2{2Ef \2mNj 
Eh {m^ - Eh 
2 (27r)^ 


/ 7Timax\2 / 7-imin\^ 

[E^ ) - [E^ ) 


(6.5) 

( 6 . 6 ) 
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Again, in the second line we have nsed the massless limit, 
rest frame can be written as 

ml + ml- ml^ 

1 

2m^ 

with similar dehnitions and kinematic bonnds, 


The nentrino energy Ei, in y’s 


(6.7) 



- 'im.y.Ek, 


now here E'^ and Ef^ are the energies of 0 and h in the rest frame, respectively. 
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